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摘 要 : 蒙古 国 风蚀 沙漠 化 严重 ,成 为 世界 主要 的 沙 尘 源 地 之 一 ,为 了 深入 了 解 蒙古 国 全 域 风 蚀 时 空 变 化 特征 , 利 


古国 境内 的 GLDAS、OLM MOD13Q1 , NASA-SRTM 等 气象 和 遥感 数据 ,基于 CEE(Google Earth Engine) 云 计算 
,利用 RWEQ 模型 计算 了 蒙古 国土 壤 风 蚀 量 和 时 空 变 化 规律 ,分析 了 土壤 风蚀 的 主要 影响 因素 及 作用 机 制 。 
果 表 明 :(1) 自 2001 一 2020 年 ,蒙古 国土 壤 风蚀 呈 北 部 低 .南部 高 的 空间 分 布 格局 ;土壤 风蚀 主要 发 生 在 春季 , 风 
蚀 量 占 全 年 的 45%; 近 20a 土 壤 风蚀 量 总 体 上 呈 显 著 增 加 趋势 ,风蚀 模 数 以 0.06 thm- a 的 速度 增长 。(2) 蒙古 国 
HE 


风蚀 变化 整体 上 与 气候 因子 及 由 其 引起 的 植被 盖 度 变化 有 关 ,与 放牧 EMIT ER OAT 3 RA SU s V] 


系 。 研 究 结果 可 为 蒙古 国 风蚀 沙漠 化 防治 提供 重要 的 科学 参考 。 


关键 词 : RWEQ 模型 ;土壤 风蚀 ;影响 因素 ; 蒙古 国 


土壤 风蚀 是 土地 退化 最 主要 的 原因 之 一 。 世 
界 上 约 三 分 之 一 的 土地 受到 风蚀 的 影响 ,其 中 严重 
风蚀 土地 占 比 达 50% 以 上 5 。 土 壤 风 蚀 会 对 陆地 
生态 系统 以 及 人 类 社会 发 展 产 生 诸多 不 利 影响 "， 
如 显著 降低 土壤 肥力 和 植物 生产 力 外 ,导致 空气 质 
量 降低 等 和 “。 此 外 ,地 表 风 蚀 形 成 的 细 粒 物质 是 大 
气 气 洲 胶 的 主要 来 源 ”, 为 沙尘暴 提供 了 物 源 。 只 
要 有 发 生 风 蚀 的 条 件 ,就 可 能 产生 较 大 规模 的 风沙 
灾害 所 ,如 2021 年 3 月 14 日 蒙古 国 发 生 的 特大 沙 尘 
暴 事件 就 危及 到 东亚 广大 地 区 。 

土壤 风蚀 的 定量 评估 方法 包括 野外 风蚀 直接 
测量 ” 、 风 洞 模拟 实验 "" .基于 同位 素 示 踪 法 的 风蚀 
速率 计算 等 看。 利用 风蚀 监测 设备 或 风 洞 模拟 实 
验 具 有 一 定 的 局 限 性 "中 ,难以 对 较 大 空间 范围 进行 
风蚀 评估 。 随 着 遥感 技术 的 发 展 , 越 来 越 多 的 学 者 
利用 地 理 信息 系统 进行 风蚀 建 模 ,发 展 出 WEQ、 
RWEQ 、WEPS .SWEEP 等 众多 模型 ,其 中 以 RWEQ 
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模型 应 用 最 为 广泛 。RWEQ 模型 以 WEOQO 模型 为 基 
础 ,融合 了 自然 过 程 和 实际 经 验 , 能 够 将 野外 数据 
与 计算 机 模型 相 结合 ,来 描述 物理 的 风蚀 过 程 拉 ， 
在 实际 中 得 到 了 广泛 的 应 用 。RWEQ 模型 起 初 是 
美国 农业 部 (USDA ) 为 估算 耕地 十 壤 流失 而 开发 的 
一 种 基于 经 验 的 模型 “ ,之 后 经 过 数 次 改进 和 验 
证 ,在 其 他 各 地 得 到 了 广泛 的 应 用 。Zhang °° 
用 RWEQ 模 型 ,对 中 国内 蒙古 地 区 1990 一 2015 年 的 
土壤 风蚀 进行 了 模拟 ,评估 了 和 气候 变化 对 风蚀 的 影 
响 ; Lin 等 52 利用 RWEQ 模型 估算 了 1982—2015 年 
中 国 河西 地 区 的 潜在 风蚀 率 ,Chi 等 %% 利 用 RWEQ 
模型 评估 了 中 国土 壤 风 人 蚀 的 时 空格 局 ,Li 等 "利用 
RWEQ 模型 模拟 了 1986—2099 年 中 亚 土壤 风蚀 的 
时 空 动态 。 蒙 古国 是 东亚 地 区 沙尘暴 事件 发 生 频 
率 最 高 的 地 区 ,是 东亚 地 区 沙尘暴 灾害 的 主要 起 
源 地 "” 。 前 人 使 用 基于 专家 的 方法 评估 了 蒙古 国 
的 风蚀 速率 和 程度 ,结果 表明 蒙古 国 大 约 30% 的 领 


Ez 


P 莹 草场 荒漠 化 防治 技术 合作 研究 与 示范 "(2017YFE0109200); "P ESSE 


P 国 科学 院 王 宽 诚 率先 人 才 计 划 “ 产 研 人 才 扶 持 项 目 " 共 同 资助 


作者 简介 : 郭 黄 (1994- ) , 女 ,硕士 研究 生 ,主要 研究 方向 为 水 土 保持 与 荒漠 化 防治 研究 . E-mail: guoyin19@mails.ucas.ac.cn 


通讯 作者 : 范 敬 龙 . E-mail: fanjl@ms.xjb.ac.cn 


http: //azr.xjegi.com 


202210.00013v1 


chinaXiv 


ChinaXiv 合 作 期 刊 
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土 受 到 风蚀 的 影响 ,其 中 戈壁 沙漠 被 认为 是 受 风 伺 
影响 最 严重 的 区 域 拉 。Natsagdorj 等 2 利用 气象 观 
测 和 激光 雷达 ,发现 蒙古 国 各 地 的 风蚀 总 体 上 有 所 
增加 。 师 华 定 等 汪 利 用 地 理 信 息 系统 和 遥感 技术 ， 
通过 模糊 ec 均 值 聚 类 (FCM ) 方 法 对 植被 覆盖 率 .地 
形 起 伏 度 .土壤 干燥 度 和 风 场 强度 4 个 主要 蒙古 高 
原 风蚀 环境 因子 数据 进行 模糊 聚 类 ,结合 专家 知 
识 , 将 环境 组 合 与 特定 的 风蚀 危险 度 相 对 应 ,得 到 
风蚀 危险 度 图 ,以 此 对 蒙古 高 原 土壤 风蚀 危险 度 进 
行 评价 研究 ,并 分 析 了 蒙古 高 原 土 壤 风 蚀 危险 度 的 
空间 格局 特征 。FCM 方 法 可 有 效 地 获取 风蚀 -环境 
关系 知识 ,为 预测 性 风蚀 危险 度 制 图 提供 依据 ,此 
研究 并 没有 对 蒙古 高 原 土壤 风蚀 危险 度 进行 定量 
评价 。Mandakh 等 3 在 ArcGIS 环境 中 使 用 风蚀 方 
程 模型 (WEQ) 评 估 了 蒙古 国 各 地 的 风蚀 情况 ,尽管 
此 研究 对 蒙古 国 的 风蚀 进行 了 定量 分 析 , 但 是 缺乏 
对 蒙古 国 长 时 间 序 列 的 风蚀 趋势 的 分 析 。 目 前 , 关 
于 蒙古 国 的 风蚀 研究 较为 有 限 , 焉 竺 开展 蒙古 国 长 
期 风蚀 状况 定量 描述 和 评估 研究。 本 文 使 用 
RWEQ 模型 ,计算 蒙古 国 全 域 的 风蚀 模 数 ,揭示 风 
乌 的 时 空 变 化 规律 和 趋势 ,分 析 识 别 主 要 影响 
素 , 有 助 于 深入 了 解 蒙古 国土 壤 风蚀 整体 特征 ,为 重 
点 地 区 风沙 灾害 和 荒漠 化 防治 工作 提供 重要 参考 。 


1 研究 区 概况 


蒙古 国 位 于 亚洲 东 中 部 (41"32' 252? 15' N, 
87°44'~119°56'E) ,面积 156.65 x10'km252 ,是 世界 
第 二 大 内 陆 国 家 。 蒙 古国 平均 海拔 约 1$80 m ,最 高 
点 (西部 的 khuiten 峰 ) 海 拔 4653 m ,最 低 点 (东部 的 
Khookh Nuur WEH ) f$ 4 518 m, AARS ŽE FET 
1542 1310 m。 蒙 古国 远离 海洋 ,具有 显著 的 大 陆 
性 气候 特征 ,四 季 分 明 , 年 均 温和 日 均 温 变化 较 大 ; 
年 平均 降水 量 约 120~250 mm, 降 水 量 从 北 到 南 .从 
东 到 西 递减 ,降水 多 集中 在 7 一 8 月 3。 根据 蒙古 国 
家 气象 和 环境 监测 局 (https://eic.mn/) 提 供 的 蒙古 国 
自然 地 带 分 布 数 据 显 示 , 南 部 为 沙漠 和 荒漠 草原 地 
区 ,东南 部 是 戈壁 沙漠 区 ;典型 草原 主要 分 布 在 莹 
古国 中 部 东西 延伸 ;和 森林 集中 分 布 在 蒙古 国 北部 山 
区 (图 1)。 其 中 ,荒漠 草原 、 半 充 漠 和 沙漠 地 区 约 占 
41.3% ,这 些 地 区 大 风 频 发 ,土壤 风蚀 发 生 强烈 ,生态 
环境 脆弱 ,年 平均 风速 为 4~6 m-s 7, Zt Edna 


0 250km 


注 : 底 图 采用 自然 资源 部 标准 地 图 制作 , 审 图 号 为 GS(2020)4395 
号 ,对 底 图 边界 无 修改 。 下 同 。 

图 1 研究 区 概况 示意 图 

Fig. 1 Geographical distribution map of Mongolia 


FS Fre oe DX p F EES ,更 易 被 风蚀 并 形成 沙 尘 
暴 ”” ,年 沙尘暴 日 数 可 达 30~120 d; 蒙 古国 是 一 个 
沙漠 化 严重 的 国家 ,沙漠 化 的 加 剧 直 接 增加 了 风蚀 
和 沙 尘 排放 源 ,威胁 到 东亚 其 他 国家 的 空气 质量 ”。 


2 研究 方法 与 数据 来 源 


2.1 研究 方法 

2.1.1 风蚀 模型 RWEQ 模型 充分 考虑 了 包括 气候 
条 件 .土壤 特性 .表面 粗糙 度 .植被 覆盖 和 其 他 环境 
要 素 ,主要 基于 式 (1)~(3) 来 估算 风蚀 量 。 


SWEQ = 2x9, ot) (1) 
S 


0.3711 


S=150.71(WF x EF x SCF x K' x COG) (2) 

Q „œx = 109.8(WF XEF XSCFXK'xCOG) (3) 

式 中 :SEO 表示 风蚀 模 数 (thm ”aa i Qua BEAN EK 

大 输 沙 能 力 (kgm');s 表示 关键 地 块 长 度 (m);x 是 模 

拟 区 域 不 可 侵蚀 边界 到 下 风口 的 距离 (m); WF 表示 

ABD T (kg ^ m); EF 表示 土壤 可 刨 性 因子 (无 量 

28) ;SCF 表示 土壤 结 皮 因子 (无 量 纲 );K' 表 示 土 壤 

粗 烟 度 影响 因子 (无 量 纲 );C0G 表 示 植 被 影响 因 
子 。 具 体 计算 过 程 如 下 式 所 示 : 


WE =W, x7 SWxSD (4) 
_W 
W,= Nd (5) 
a 2 
W- Y UU, -U)) (6) 
i=l 


s 348.0 1013 -0.1 E 1.004882" (7) 


SD = 1 - P(snow cover > 25.4 mm) (8) 
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式 中 :和 % 是 第 ;年 和 第 /年 的 风蚀 模 数 值 或 风速 、 


E (kg: m’); g 表示 重力 加 速度 (ms ;SW 表示 土壤 


湿度 因子 ;SD 表示 积 雪 覆盖 影响 因子 ; 玉 表 示 风 力 
[E (m? - 7) ; Usd 2 m 高 的 风速 (ms ;以 指 风速 效 
值 ,一 般 取 5 ms';Nd 是 风速 测量 的 时 间 间 隔 (d); 
NN 代表 风速 观测 频率 (d) ;EL 表示 从 数字 高 程 模型 
(DEM ) 数 据 获得 的 高 程 (km) ;7 表示 绝对 温度 (K); 
SD 由 积 雪 和 宪 盖 度 计 算 获 得 ;P 表 示 测 定时 段 内 积 雪 
履 盖 深度 大 于 25.4 mm 的 发 生 概 率 。 
土壤 可 蚀 性 因子 EF YET SERE PAIR : 
EF =(29.09 + 0.31Sa  0.17Si + 0.3379. 2 
2.590M - 0.95CaCO,)/100 
式 中 :Sa 为 土壤 砂粒 含量 (%);S 为 土壤 粉 粒 含量 
(96); Cl lh E t (96) 5 Sal CI H HED BS b b 
Lk; OM 为 土壤 有 机 质 含 量 (%o); CaCO: 为 碳酸 钙 含量 
(9%)。 
土壤 结 皮影 响 因 子 SCP 计算 过 程 如 下 式 所 示 : 


SCF= 1 = (10) 
1 +0.0066(Cly + 0.021(0MY 


土壤 质地 和 有 机 质 含量 的 年 际 变化 非常 缓慢 ， 
可 将 土壤 可 蚀 性 因子 和 结 皮 因 子 视 为 常量 值 , 假 设 
土壤 可 蚀 性 因子 和 结 皮 因 子 不 随时 间 改 变 。 
地 表 粗 糙 度 影响 因子 及 ,根据 Xu 等 后 的 研究 
结果 ,计算 如 下 式 所 示 : 
K'=cosa (11) 
式 中 :a 表 示 地 形 坡 度 , 可 从 DEM 数据 中 提取 。 
植被 影响 因子 COC 表 示 一 定植 被 覆盖 条 件 下 
对 土壤 风蚀 的 抑制 程度 ,采用 最 大 合成 法 (MVC ) 得 
到 半月 NDVI 数 据 , 并 用 像 元 二 分 法 求 得 2001 一 
2020 年 的 半月 植被 覆 善 度 ,计算 过 程 如 下 式 所 示 : 
COC=e 0 (12) 
NDVI- NDVI 
NDVI,„ -NDVI 
式 中 :FVC 为 植被 覆盖 率 (%);NDVI 为 宰 土 像素 的 
NDVI 值 ;NDVI 为 植被 的 NDVI 值 。 
2.1.2 交 化 趋势 分 析 方 法 ”利用 Sen’s 趋 热 分析 法 ， 
计算 了 2001 一 2020 年 蒙古 国土 壤 风蚀 以 及 气象 因 
子 (风速 气温、 降水 ) 的 变化 趋势 。 该 方法 是 非 参 
数 检验 ,只 需要 保证 样本 的 独立 性 ,不 需要 拟 合 正 
dfi HR PAIR: 


(9) 


FVC- (13) 


(s -x) 
Slope = Median E (Vi «isj&n) (14) 


l 


气温 、 降 水 的 数值 ;n 是 研究 期 的 总 年 数 。Slope 表 示 
时 间 序 列 数据 的 变化 率 , 正 值 表示 风蚀 或 气象 因子 
(风速 气温、 降水 ) 随 时间 变化 总 体 呈 增加 趋势 , 负 
值 表示 总 体 呈 减少 趋势 。 
2.2 数据 来 源 

(1) 气象 数据 

气象 数据 包括 风速 .气温 、 降 水 、 积 雪 深 度 、 土 
壤 水 分 , 均 来 自 全 球 陆地 数据 同化 系统 GLDAS2.1 

品 。GLDAS 数 据 是 美国 航空 航天 局 (NASA ) 戈 

达 德 空间 飞行 中 心 (GSEC) 和 美国 海洋 和 大 气 局 
(NOAA) 国 家 环境 预报 中 心 (NCEP) 联 合 发 布 的 基 
于 卫星 . 陆 面 模式 和 地 面 观测 数据 的 同化 产品 , 它 
能 够 提供 多 种 驱动 数据 ,这 些 数据 来 源 于 大 气 同化 
产品 .再 分 析 和 实际 观测 。GLDAS2.1 产 品 时 间 分 辩 
KH 3 h, ZS [8] 4) 9t 9 2j 0.25? ,与 地 面 气象 观测 站 实 
测 数据 具有 较 高 的 一 致 性 ?”"。 

(2) 土壤 数据 

土壤 数据 包括 土壤 有 机 质 含 量 . 黏 粒 含量 、 粉 
粒 含量 、 沙 粒 含量 .CaC0; 含 量 , 其 中 土壤 有 机 质 含 
AE a 、 粉 粒 含 量 、 沙 粒 含量 来 源 于 Open 
Land Map (OLM) 数 据 集 ,CaC0; 含 量 来 源 于 世界 土壤 
数据 库 (HWSD )。 

(3) 植被 数据 

选用 由 美国 航空 航天 局 (NASA ) 发 布 的 NDVI 
中 分 状 率 成 像 光谱 仪 获得 的 陆地 2 级 标准 数据 产品 


MOD13Q1, 时 间 分 辩 率 为 16 d, 空间 分 辨 率 为 250 
m, 可 用 于 计算 植被 覆盖 度 。 
(4) 地 形 数 据 


选择 由 美国 航空 航天 局 (NASA) 和 美国 测绘 局 
(CNIMA) 联 合 测量 的 数字 高 程 模 型 (SRTM) 版 本 4， 
空间 分 辩 率 为 90 m。 这 是 目前 质量 最 高 的 SRTM 数 
据 集 ,水 平 标准 误差 为 1 m, 垂 直 标 准 误差 为 0m ,可 
用 于 计算 坡度 和 地 表 粗 糙 度 。 

(5) 其 他 数据 

牲畜 数量 数据 来 自 于 蒙古 国 国 家 统计 局 
(www.1212.mn)。 土 地 利用 数据 使 用 的 是 由 中 
家 基础 地 理 信息 中 心 (NGCC) 提 供 的 GlobeLand30 
产品 中 ,数据 集 使 用 了 包括 20000 多 张 美国 陆地 资 
源 卫 星 (Landsat) 和 中 国 环境 减灾 卫星 (HJ-1) 的 多 
光谱 影像 ,具有 和 较 高 的 空间 分 辨 率 和 数据 精度 。 

将 以 上 所 有 栅 格 数据 的 空间 分 辨 率 重 新 采样 
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至 500 m, 时 间 分 状 率 为 16 d, 全 年 有 23 个 时 间 段 。 
对 计算 结果 进行 汇总 ,得 到 每 年 的 风蚀 总 量 。 


3 BRAD 


3.1 风蚀 模 数 的 空间 变化 特征 

从 蒙古 国土 壤 风 人 蚀 模 数 空间 分 布 来 看 (图 2)， 
土壤 风蚀 变化 空间 差异 较 大 ,风蚀 模 数 的 最 大 值 与 
最 小 值 的 差 值 在 0~60 thm2.a-: 之 间 。 整 体 上 , 土 
坏 风 蚀 模 数 呈现 北部 低 .南部 高 的 空间 分 布 模式 ， 
由 于 南部 的 戈壁 沙漠 区 土壤 质地 较 粗 , 物 源 丰富 ， 
且 可 蚀 性 较 高 ,加 之 降水 少 ,地 表 植 被 稀 玻 ,因此 该 
区 成 为 蒙古 国土 壤 风蚀 最 为 严重 的 地 区 。 


(a) 2001 年 (b) 2002 年 


SWEQ(t-hm?-a") 


(c) 20034 


由 图 3 .图 4 可 知 , 对 风蚀 模 数 进 行 多 年 平均 ,从 
年 内 空间 分 布 可 以 看 出 ,风蚀 整体 上 春季 (3 一 5 月 ) 
最 强 , 约 占 全 年 风蚀 模 数 的 45% ; 其 次 为 秋季 (9 一 
11 月 ) , 约 占 全 年 风蚀 模 数 的 33% ;再 次 为 夏季 (6 一 
8 月 ) , 约 占 全 年 风蚀 模 数 的 12% ;冬季 (12 月 一 次 年 
2 月 ) 最 弱 , 约 占 全 年 风蚀 模 数 的 10%。 和 上 蒙古 高 原 春 
季 降 水 稀少 ,土壤 干燥 ,植物 尚未 进入 快速 生长 期 ， 
植被 覆盖 度 低 ,因此 春季 风蚀 非常 强烈 。 冬 季 土 壤 
表层 冻结 , 且 多 有 积 雪 和 覆盖 ,因此 风 钟 最 弱 。 夏 季 
和 秋季 植被 条 件 均 较 好 ,但 夏季 降水 集中 ,土壤 湿度 
明显 大 于 秋季 ,因此 夏季 风蚀 强度 小 于 秋季 。 和 气象 观 
测 表 明 ,蒙古 国 的 沙 尘 事件 主要 发 生 在 春季 5 ,与 本 
文风 蚀 主要 发 生 在 春季 的 结论 一 致 。 


(d) 2004 年 


图 2 2001 一 2020 年 蒙古 国土 壤 风蚀 模 数 空间 分 布 特征 


Fig.2 Spatial distribution characteristics map of wind erosion in Mongolia from 2001 to 2020 
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图 3 2001 一 2020 年 蒙古 国 季 节 平 均 土 壤 风 蚀 模 数 空间 分 布 


Fig.3 Spatial distribution map of seasonal average soil wind erosion in Mongolia from 2001 to 2020 


JA piu C hm:?-a71) 


0 
2002 2004 2006 2008 2010 2012 2014 2016 2018 2020 
年 份 


图 4 2001 一 2020 年 蒙古 国 季节 平均 土壤 风蚀 模 数 变化 
Fig.4 Variation of seasonal average soil wind erosion 
modulus in Mongolia from 2001 to 2020 


3.2. 风蚀 模 数 的 时 间 变 化 特征 

由 图 5 可 知 ,在 2001 一 2020 年 间 ,蒙古 国土 壤 
风蚀 模 数 总 体 呈 显著 上 升 趋势 (已 =0.2258, P< 
0.05) ,其 中 2018 年 最 强 (4.4463 thm? a^) ,2008 年 
15:95 (1.369 thm 2.a )。 从 时 间 序 列 的 变化 特征 来 
看 ,2005 一 2009 年 的 风蚀 区 域 明 显 小 于 其 他 年 份 ， 
而 2009 年 后 风蚀 范围 明显 增 大 。 根 据 风 蚀 模 数 的 
年 际 变化 ,2001 一 2020 年 期 间 风 蚀 变 化 分 为 2 个 明 
显 的 阶段 ,2001 一 2008 年 为 风蚀 模 数 下 降 时 期 (R= 
0.536, P«0.05) , 2010—2020 年 为 风蚀 模 数 增加 
时 期 。 

从 土壤 风蚀 变化 速率 空间 差异 来 看 (图 6) , 蒙 


古国 南部 戈壁 沙漠 地 区 的 风蚀 变化 速率 约 为 21. 
hm? ar! ,明显 高 于 其 他 区 域 ;蒙古 国 北 部 风蚀 最 低 ， 
变化 最 小 ;蒙古 国 西 南部 风蚀 速率 变化 为 负 , 表 明 
此 区 域 风蚀 呈 下 降 趋 势 ,平均 每 年 的 风蚀 变化 速率 
约 -0.01~-0.6t:hm ”a"。 总 体 来 看 ,蒙古 国土 壤 风 
蚀 模 数 高 值 区 变化 速率 快 ,而 低 值 区 变化 速率 慢 。 
蒙古 国 北 部 以 森林 生态 系统 为 主 ,植被 类 型 以 高 山 
泰 加 林 高山 森 林 和 森林 草原 为 主 , 植 被 覆盖 度 较 
高 ,土壤 风蚀 阔 值 不 大 ,整体 森林 生态 系统 稳定 , 受 
气候 变化 及 人 为 活动 影响 小 ,风蚀 很 小 , 且 变 化 率 
很 小 ;蒙古 国 南 部 荒漠 生态 系统 脆弱 ,植被 稀少 , 受 
气候 变化 影响 大 ,因此 风蚀 年 际 变化 率 也 较 大 ; 蒙 
十 国 西南 部 风蚀 呈 下 降 趋势 的 主要 原因 是 此 区 域 
风速 星 下 降 趋 势 , 风 蚀 量 常年 较 低 。 

3.3 土壤 风蚀 的 影响 因子 分 析 

3.3.1 气候 因子 变化 及 影响 “影响 土壤 风蚀 的 主要 
气候 因子 有 风速 气温、 降水 等 。 本 文采 用 Sen’s 趋 
势 分 析 法 分 析 了 蒙古 国 主要 气候 因子 变化 趋势 及 
空间 分 布 。 从 图 7 可 以 看 出 , 近 20 a 来 蒙古 国 年 最 
大 风速 呈 明 显 波 动 变 化 ,2008 年 最 低 (14.04 m-s”), 
2001 年 最 高 (16.26 m * s) ; 2001—2008 4E. , 4E X 
PEEL FRAGE S (920.427 , P0.05) ; 2009—2020 4E , 4E 
最 大 风速 旦 显著 增加 趋势 (R=0.3959,P<0.05)。 土 
壤 风 蚀 模 数 随 着 风速 变化 趋势 的 波动 而 波动 ,二 者 
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45 | 一 风蚀 模 数 
| —- 2001 一 2008 年 变化 趋势 
so. = 2009_2020 年 变化 趋势 y=0.0624x-122.1049 
pi — 2001 一 2020 年 变化 趋势 R?-0.2258, P«0.05, n=20 
W 3.5 H 
D 上 
= 30 Fed Sb LLL ZO 
Bos PA 5 NN ur, 
B amm 
K 2.0 上 -一 d y=0.0905x-179.5113 
us | 25 70.1773x+357.6574 in d R’-0.194, P20.05, n-12 
7 | R-0.536, P<0.05, n=8 
S8 82 S333332 8283333333333233 
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图 5 2001 一 2020 年 蒙古 国 多 年 平均 土壤 风蚀 模 数 变化 趋势 
Fig.5 Variation trend graph of average soil wind erosion in Mongolia from 2001 to 2020 


风蚀 变化 速率 /thm?.a0 
mm 一 0.8~ -0.1 = 0.1-0.5 
mu —0.1- —0.01 mm 0.5~1 
E3-001-0 m 1-27 
= 0-0.1 


图 6 2001 一 2020 年 蒙古 国土 壤 风蚀 模 数 变化 趋势 的 空间 分 布 


Fig. 6 Spatial distribution of variation trend of soil wind erosion in Mongolia from 2001 to 2020 
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图 7 2001 一 2020 年 蒙古 国 年 最 大 风速 变化 趋势 


Fig. 7 Variation trend of annual average maximum wind speed in Mongolia from 2001 to 2020 
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变化 规律 基本 呈现 对 应 关系 (图 5, 图 7) ,表明 风速 
是 影响 土壤 风蚀 最 重要 的 自然 动力 条 件 。 

结合 年 均 风 速 变化 趋势 空间 分 布 (图 8) 来 看 ， 
蒙古 国 南部 戈壁 沙漠 区 风速 显著 上 升 , 其 风蚀 量 呈 
增加 趋势 ;中 部 风速 下 降 ,风蚀 量 呈 减 少 趋势 。 此 
外 ,将 蒙古 国 2001 一 2020 年 逐年 平均 风速 与 逐年 风 
蚀 模 数 进行 相关 性 分 析 ( 图 9) 发 现 , 在 蒙古 国 仅 
3.3% 左 右 的 区 域 平均 风速 与 风蚀 模 数 旺 负 相 关 关 
系 , 绝 大 数 区 域 (96.6% ) 风 蚀 模 数 都 与 平均 风速 呈 
正 相 关 。 尤 其 是 在 蒙古 国 南部 及 西南 部 的 戈壁 区 
域 ,风速 与 风蚀 量 呈 显著 正 相关 关系 ,说 明 风 速 是 
影响 土壤 风蚀 最 重要 的 自然 动力 条 件 。 这 也 解释 
了 蒙古 国 西南 部 风蚀 量变 化 速率 下 降 是 受 风速 下 
降 的 影响 。 

气温 和 降水 会 通过 影响 土壤 湿度 和 植物 生长 
而 间接 影响 风蚀 3。 从 图 10 可 以 看 出 , 近 20 a 来 ， 


图 例 


cx 湖泊 风速 ms ) 
L^ m. 
0 250km 


图 8 2001 一 2020 年 蒙古 国 年 均 风 速 变化 趋势 空间 分 布 
Fig. 8 Spatial distribution of annual average wind speed 
change trend in Mongolia from 2001 to 2020 


蒙古 国 大 部 分 地 区 年 均 温 呈 上 升 趋势 , 仪 有 北部 少 
部 分 区 域 年 均 温 呈 下 降 趋 势 。 南 部 增 温 最 高 , 增 温 
速率 达到 0.06 C «a^ ,升温 促进 了 地 表 水 的 蒸发 , 导 
致 土壤 湿度 降低 ,可 蚀 性 增加 ,特别 是 在 缺乏 植被 
的 地 区 ,升温 增加 了 风蚀 的 发 生 。 


相关 系数 
m —0.4~ -0.2 
mm -0.2~0 
= 0-04 

= 0.4-0.7 
mm 0.7-0.96 


图 9 蒙古 国 年 均 风 速 与 风蚀 模 数 相关 关系 


Fig. 9 Correlation between annual average wind speed and wind erosion modulus in Mongolia 
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图 10 2001 一 2020 年 蒙古 国 年 均 温 变化 趋势 时 空 分 布 


Fig. 10 Temporal and spatial distribution of annual mean temperature change trend in Mongolia from 2001 to 2020 
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降水 与 表层 土壤 含水 量 密切 相关 ,降水 量 减少 
增加 了 风蚀 的 发 生 。 整 体 来 看 , 近 20 a 来 蒙古 国 降 
水 呈 显 著 增 加 趋势 (图 11b) ,但 降水 变化 存在 较 大 
空间 差异 (图 11a) ,北部 地 区 降水 明显 增加 ,中 部 和 
南部 部 分 地 区 降水 减少 ,尤其 南部 戈壁 地 区 呈 明 显 
减少 趋势 (图 11c)。 南 部 地 区 干旱 灾害 发 生 的 可 能 
性 增加 ,加 上 风力 的 作用 ,风蚀 量 也 呈 增 大 趋势 。 
3.3.2 植被 履 盖 度 变 化 及 影响 “植被 覆盖 可 增加 地 
表 粗 糙 度 和 削弱 近 地 表 风速 ,从 而 降低 地 表 风 蚀 危 
PES), R12 为 近 20 a 来 蒙古 国 FVC 的 变化 ,可 
以 看 出 ,蒙古 国 大 部 分 区 域 FVC 整 体 呈 增长 趋势 ， 
东部 部 分 地 区 植被 覆盖 度 快 速 增 加 ,因此 土壤 风蚀 
降低 ;但 部 分 地 区 尤其 是 南部 和 西部 的 车 漠 草 原 区 
多 地 植被 覆盖 度 减 少 ,导致 风蚀 发 生 风 险 增 大 。 


4 讨论 


4.1 风蚀 模 数 的 结果 验证 

RWEQ 模型 起 初 是 美国 农业 部 (USDA ) 为 估算 
耕地 土壤 流失 而 开发 的 一 种 基于 经 验 的 模型 王 ,之 
后 经 过 数 次 改进 和 验证 ,不 仅 是 将 其 用 于 评估 农田 
发 生 的 土壤 风蚀 ,如 Li 等 将 其 用 于 评估 内 蒙古 地 
区 的 固沙 服务 功能 , Chi 等 9 评估 了 中 国 全 域 的 土 
壤 风 刨 ,也 有 一 些 学 者 将 其 用 于 中 亚 地 区 的 土壤 风 
刨 评 做 。 然 而 其 结果 往往 需要 其 他 方法 售 算 的 
结果 来 对 比 验 证 。 


(a) 蒙古 国 年 均 降水 量 空间 变化 


综 上 所 述 , 蒙 古国 2001 一 2020 年 多 年 平均 土壤 
风蚀 模 数 为 2.736 t-hm® a'o "Css Exi de fh 53 E 
壤 侵 蚀 速 率 的 常用 方法 。 前 人 基于 Cs 监测 方 
法 计算 了 蒙古 高 原 的 风蚀 模 数 ,蒙古 国境 内 的 监测 
样 点 从 北向 南 涵盖 了 和 森林 草原 区 .典型 草原 区 SETS 
化 草原 区 .草原 化 荒漠 区 .戈壁 区 BUS EDO 
将 本 人 研究 中 基于 RWEQ 模型 模拟 的 20 a 的 土壤 风 
蚀 模 数 与 刘 纪 远 等 ”使 用 Cs 实测 的 土壤 风蚀 模 
数 基于 相同 的 地 理 坐 标点 进行 了 相关 性 分 析 ( 图 
13) ,结果 表明 ,RWEQ 模 型 模拟 结果 与 ?Cs 测量 值 
显著 相关 (R=0.91, P«0.01) Mandakh 等 后 使 用 
WEQ 模 型 评估 的 蒙古 国 风 蚀 模 数 在 2.7~27.5 te hm: 
a zz [RI , 比 本 文 用 RWEQ 估算 结果 略 偏 低 ,但 数量 
值 大 体 上 相当 。 土 壤 风蚀 的 发 生 和 发 展 依赖 于 侵 
蚀 动 力 (气候 因子 如 风速 ) 和 可 蚀 性 因子 (地 表 状 况 
如 植被 盖 度 .土壤 特性 等 ) 之 间 的 相互 作用 ,土壤 可 
蚀 性 因子 可 以 衡量 土壤 表层 被 风力 侵蚀 的 相对 可 
能 性 。Jugder 等 中 从 土壤 可 刨 性 的 角度 评价 了 蒙古 
国土 壤 被 风力 侵 刨 的 相对 可 能 性 ,发 现 从 蒙古 国 南 
部 的 沙漠 和 荒漠 草原 地 区 到 西部 的 沙 地 ,存在 着 广 
泛 的 高 可 蚀 性 区 域 , 与 本 文中 的 严重 风蚀 区 域 也 大 
致 相同 。 中 等 可 蚀 性 区 域 主要 覆盖 西部 的 干旱 草 
原 区 和 荒漠 草原 区 , 低 可 蚀 性 区 域 从 蒙古 国 西部 延 
伸 到 北部 。Mandakh EEH WEQ 模型 评估 了 蒙 
古国 风蚀 模 数 ,分 析 了 风蚀 的 空间 分 布 规律 , 与 本 
文 结 论 一 致 。 


(6) 蒙古 国 年 均 降水 量 时 间 变 化 


y=6.0454x-11904 
R?=0.0159, P<0.01 


年 均 降 水 量 /mm 


2000 2005 2010 2015 2020 
年 份 


(c) 蒙古 国 南 部 年 均 降水 量 时 间 变 化 


y- -0.6514x-1442 
175 R= 0.0159, P>0.5 
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125 
100 
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图 11 2001 一 2020 年 蒙古 国 年 降水 量变 化 趋势 时 空 分 布 


Fig. 11 Temporal and spatial distribution of annual rainfall variation trend in Mongolia from 2001 to 2020 
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FVC 变 化 


口 省 界线 = -0.03~ —0.001 
= —0.001-0 

P en = 0-0.002 
= 0.002-0.01 
= 0.01~0.041 


图 12 2001 一 2020 年 蒙古 国 植被 覆盖 度 变 化 趋势 


Fig. 12 Change trend of annual vegetative cover in Mongolia 
from 2001 to 2020 
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13 ”Cs 测量 和 RWEQ 模 型 模拟 风蚀 模 数 的 相关 系数 
Fig. 13 Correlation coefficient of ”Cs measurement and 


RWEQ model simulation wind erosion modulus 


42. 人 为 活动 因素 变化 及 影响 
土壤 风蚀 的 发 生 与 气候 条 件 和 下 垫 面条 件 有 


图 例 
mu 草地 转 耕 地 
国 草地 转 林 地 
mu 草地 转 建设 用 地 
mu 未 利用 地 转 草地 MN 林地 转 裸 地 
mu 未 利用 地 转 耕 地 EE 林地 转 草地 
mm 未 利用 地 转 林地 ES 林地 转 耕 地 


密切 关系 ,人 为 活动 通过 影响 下 扑面 而 影响 风蚀 的 
发 生 。 过 度 放牧 以 及 煤矿 开采 等 人 类 活动 加 剧 了 
蒙古 国土 地 退化 下 。 由 图 14 可 知 , 在 2001 一 2020 
年 间 , 蒙 古国 的 牲畜 数量 由 2001 年 的 26075 千 头 增 
加 到 2020 年 的 67068 千 头 , 每 年 大 约 以 2400 千 头 的 
速度 显著 增加 。 牲 冀 数 量 的 增加 使 草原 压力 增 大 ， 

宽松 的 放牧 政策 造成 多 数 草场 过 度 放 牧 ,牲畜 踩 
踏 和 过 牧 行为 破坏 了 土壤 表层 结构 ,降低 了 冬 E 
季 地 表 的 植被 盖 度 ,增加 了 土壤 的 风蚀 风险 。 
此 ,过 度 放牧 可 能 是 导致 蒙古 国土 壤 风 刨 增加 的 一 
个 人 为 活动 因素 。 
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图 14 2001 一 2020 年 蒙古 国 牲畜 数量 变化 趋势 
Fig. 14 Variation trend of livestock quantity in Mongolia 
from 2001 to 2020 


UE 20 a 来 ,蒙古 国 经 济 社会 快速 发 展 ,为 了 适应 
生存 和 发 展 的 需求 ,土地 利用 类 型 也 发 生 了 较 大 变 
化 。 从 土地 利用 类 型 转化 (图 15 ) 来 看 ,蒙古 国境 内 


mu 耕地 转 草地 

国 草地 转 未 利用 地 

mm 未 利用 地 转 建设 用 地 
um 其 他 

C1 未 变化 


图 15 2001 一 2020 年 蒙古 国 主要 土地 利用 类 型 转化 
Fig. 15 Transformation of main land use types in Mongolia from 2001 to 2020 
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有 1733 km 的 草地 转化 为 耕地 ,部 分 未 利用 地 和 林 
地 也 有 转化 为 耕地 的 情况 。 与 传统 畜牧 业 生 产 方 
式 相 比 ,发 展 种 植 业 ,会 导致 更 严重 的 土壤 风蚀 发 
生 。 其 他 土地 类 型 转化 为 耕地 后 ,人 为 翻动 表土 ， 
土壤 结构 水 分 和 养分 状况 发 生 了 很 大 变化 ,导致 
局 部 地 区 风蚀 加 剧 。 
4.3 影响 因素 耦合 的 空间 异 质 性 作用 

气温 和 降水 通过 影响 土壤 湿度 和 植被 覆盖 度 
而 间接 地 影响 风蚀 。2001 一 2020 年 蒙古 国 全 域 气 
温和 降水 均 呈 增加 趋势 ,植被 覆盖 度 也 呈 增 加 趋 
势 , 对 地 表 风 蚀 应 发 挥 一 定 的 抑制 作用 ,但 全 域 风 
蚀 仍 整体 呈 增 加 趋势 ,似乎 存在 一 定 的 矛盾 。 但 分 
析 发 现 , 原 因 在 于 影响 风蚀 因素 耦合 存在 较 大 的 时 
空 异 质 性 。 对 于 一 些 风蚀 相对 严重 的 地 区 ,植被 覆 
盖 度 的 增加 对 减少 风蚀 没有 明显 影响 ,这 是 因为 高 
风蚀 区 域 的 植被 覆盖 度 较 低 , 上 且 增 加 不 显著 (图 12)， 
风速 等 其 他 因素 的 变化 对 风蚀 有 更 显著 的 影响 ,而 
且 在 正 值 风 季 的 春季 作用 较为 强烈 。 学 者 们 对 蒙古 
国土 壤 风蚀 的 研究 中 发 现 ,不 同 的 气候 因素 影响 不 
同 的 区 域 ,干旱 对 戈壁 和 沙漠 地 区 影响 更 大 咏 。 
此 ,从 蒙古 国 全 域 的 尺度 来 看 ,风速 是 土壤 风蚀 发 
生 的 主导 因素 ,风速 越 大 ,风蚀 越 强 。 而 且 强 风蚀 
的 干旱 荒漠 土地 面积 占 比 极 高 ,因此 蒙古 国 整体 风 
蚀 仍 呈 增 加 趋势 。 


5 结论 


本 文通 过 RWEQ 模 型 对 蒙古 国 2001 一 2020 年 
的 土壤 风蚀 时 空 变化 状况 及 对 主要 影响 因子 进行 
了 分 析 , 得 到 的 结论 如 下 : 

(1) 蒙古 国土 壤 风 蚀 年 内 季节 变化 明显 ,全 年 
45% 的 风蚀 发 生 在 春季 ,秋季 其 次 ,冬季 最 少 。 

(2) 从 2001 一 2020 年 ,蒙古 国土 壤 风 蚀 呈 年 际 
波动 变化 ,总 体 上 呈 显 闭 增 加 趋势 ,风蚀 模 数 以 
0.06 t-hm?-a APR HEX s 

(3) 土壤 风蚀 具有 明显 的 空间 异 质 性 ,整体 呈 
南部 高 .北部 低 的 空间 分 布 格局 ,风蚀 严重 区 域 位 
于 蒙古 国 南部 与 中 国 接壤 的 戈壁 和 沙漠 地 区 ;风蚀 
变化 速率 空间 变异 较 大 ,土壤 风蚀 高 值 区 大 于 低 
值 区 。 

(4) 蒙古 国土 壤 风蚀 变化 整体 上 与 气候 因子 及 
引起 的 植被 盖 度 变化 有 关 ,也 与 过 牧 、 耕 地 开县 等 
人 为 活动 增强 有 密切 关系 。 为 了 有 效 防治 沙漠 化 


危害 ,蒙古 国 应 在 土地 资源 开发 利用 方面 加 强 科 学 
管理 ,对 于 严重 沙漠 化 地 区 应 积极 采取 可 行 的 治理 
措施 。 
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Soil wind erosion characteristics and main influencing factors in Mongolia 


in recent 20 years 
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Abstract: The accurate evaluation of wind erosion dynamics is important to reduce the damage caused by wind 
erosion in Mongolia. In this study, we used GLDAS, OLM, MOD13Q1, NASA-SRTM, and other meteorological 
and remote sensing data and the revised wind erosion equation model on the Google Earth Engine cloud 
computing platform to explore the spatio-temporal dynamics of soil wind erosion in Mongolia from 2001 to 2020. 
Results showed that wind erosion modulus had a significantly increased trend during the study period. In general, 
soil wind erosion tended to increase significantly over the last 20 years with the modulus of wind erosion 
increasing at a rate of 0.06 t- hm .a '. Spatially, wind erosion was strong in the southern regions and weak in the 
northern regions. Meanwhile, wind erosion in spring accounted for 4596 of the whole year. Changes in soil wind 
erosion in Mongolia were by and large identified with climatic factors and changes in vegetation cover, but they 
were also firmly identified with expanded anthropogenic exercises, for example, overgrazing and development of 
arable land. This study will provide a comprehensive understanding of the development trend of wind erosion 
desertification in Mongolia and a reference value to study dust storm in East Asia. 


Keywords: RWEQ model; wind erosion; influencing factors; Mongolia 


